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We study the Cooper minimum in high harmonic generation from argon atoms using long wave- 
length laser pulses. We find that the minimum in high harmonic spectra is systematically shifted 
with respect to total photoionization cross section measurements. We use a semi-classical theoreti- 
cal approach based on Classical Trajectory Monte Carlo and Quantum Electron Scattering methods 
(CTMC-QUEST) to model the experiment. Our study reveals that the shift between photoioniza- 
tion and high harmonic emission is due to several effects: the directivity of the recombining electrons 
and emitted polarization, and the shape of the recoUiding electron wavepacket. 



I. INTRODUCTION 

High order harmonic generation (HHG) takes place 
when an atom or a molecule is submitted to a strong 
laser field (with intensities / in the range 10^^ — 10^^ 
W.cm~^). Under the influence of the strong field, an elec- 
tron can be tunnel-ionized, accelerated and driven back 
to its parent ion where radiative recombination results 
in the emission of XUV (extreme ultraviolet) radiation 
[J, 0]. While this process was initially considered as a 
secondary source of radiation, it was soon realized that 
it is closely related to photoionization and that it could 
thus encode structural information on the irradiated tar- 
get. The first experimental observation that confirmed 
this link has consisted of the appearance of a local mini- 
mum in the high harmonic emission from argon atoms 3] , 
which was associated to the Cooper minimum observed 
in XUV photoionization of argon Q. In photoioniza- 
tion it is well known that this minimum is due to a zero 
dipole moment between the p ground state wavefunction 
and the d wavefunction of the photoionized electron for a 
photon energy about 48 eV. Such minima have been ex- 
tensively studied in photoionization because they consti- 
tute clearly identifiable features against which theoretical 
models can be tested. 

The interest for structural minima in HHG has been re- 
vived when the case of high order harmonics from aligned 
molecules was considered, showing that minima encoding 
the molecular structure could appear in harmonic spec- 
tra [H] ■ Many works have focused from then on observing 
these minima These studies have raised a number 

of questions on the modelling of the HHG process such 
as the influence of the ionic potential on the recolliding 
electron jl4Lll5l|. t he contribution of multiple molecular 
orbitals |lG.[l6lIl7|. the influence of the strong laser held 
[H, [l^ and the role of propagation effects [20]. Since 
some of these questions are still debated, and since HHG 
from molecules is quite complex to model, we decided to 



come back to the simpler case of high order harmonics 
from atoms. 

In this paper, we study the Cooper minimum in high 
order harmonic emission from argon atoms using tun- 
able infrared (1800-2000 nm) femtosecond laser pulses. 
We perform a systematic experimental study of the po- 
sition of the minimum as a function of the laser held 
(intensity and wavelength) and macroscopic parameters 
(gas pressure, beam focusing conditions). We flnd a sys- 
tematic shift of more that 5 eV in the position of the 
minimum with respect to total photoionization cross sec- 
tion measurements [111, We perform a theoretical 
study to understand the origin of this shift and flnd 
that it is partly due to the difference between angle- 
integrated photoionization measurements with unpolar- 
ized light and the HHG recombination process which is 
highly differential with respect to both electron and po- 
larization directions. Using a semi-classical simulation 
based on a combination of Classical Trajectory Monte 
Carlo (CTMC) [23] and Quantum Electron Scattering 
techniques (QUEST), we show that the additional con- 
tribution to the shift is due to the shape of the recolliding 
electron wavepacket. This new theoretical description, to 
which we refer to as CTMC-QUEST properly accounts 
for the influence of the ionic potential on the recolliding 
electron wavepacket and is thus able to reproduce accu- 
rately the experimental high order harmonic spectra. 

The manuscript is organized as follows. In section II, 
we describe the experimental setup and present the mea- 
surements of the Cooper minimum in various laser and 
macroscopic conditions. In section HI, we analyze the 
analogies and differences between photoionization and 
recombination using quantum electron scattering theory. 
Finally we perform in section IV a complete simulation 
of the experiment combining CTMC and QUEST ap- 
proaches. Atomic units are used throughout the paper 
unless otherwise stated. 
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FIG. 1: Experimental setup (see text). 



II. EXPERIMENT 



camera (PCO). 

An important item in the experimental setup is a 200 
nm-thick aluminum filter placed before the spectrome- 
ter. Indeed, high harmonic radiation produced by the 
1800 nm pulses can extend to more than 100 eV. Our 
study deals with measurements of the Cooper minimum, 
around 50 eV. Since the high harmonic emission from ar- 
gon is quite efRcient around 100 eV while it is minimum 
around 50 eV [131, the second order diffraction of the 
100 eV radiation by the grating can significantly affect 
the shape of the spectrum at 50 eV. We have observed 
this effect as a bump which fills in the Cooper minimum 
as the laser energy increases. The aluminum foil in our 
experiment filters out the radiation above 73 eV, enabling 
us to get rid of this artifact. 



A. Experimental setup 

It is well known [24] that harmonic spectra generally 
consist of a rapid decrease of the yield for low harmonic 
orders, followed by a plateau which extends up to the 
cut-off photon energy ~3.17?7p + Ip, where Up is the 
laser ponderomotive energy Up = I/Aloq, with ujq the 
laser frequency, and Ip is the target ionization potential 
{Ip =15.76 eV for Ar). Therefore, an accurate determi- 
nation of the position of the Cooper minimum in HHG 
can be achieved provided the minimum belongs to the 
plateau region of the spectrum. Since the minimum lies 
around 50 eV, this requires using laser intensities above 
2 X 10^^ W.cm~^, close to the saturation intensity of ar- 
gon when using 800 nm ^ 40 fs laser pulses. A way to 
perform measurements at higher intensities is to use few- 
cycle laser pulses: Worner et al. recently measured the 
Cooper minimum in Argon at / up to 3.5 x 10^"* W.cm~^ 
using 8 fs 780 nm laser pulses ^25|]. We have chosen an 
alternative solution to obtain a high cut-off: using long 
wavelength driving laser [261-5^. The cutoff frequency 
scales as / X A^, making it possible to produce 100-eV 
photons at less than 1 x lO^'' W.cm^^ using 1800-nm fs 
pulses. 

The experimental setup is presented in Fig. [T] We use 
the IkHz Aurore laser system from CELIA which delivers 
800 nm, 7 mJ, 35 fs pulses. We reduce the pulse energy 
to 4.7 mJ and inject it to a HE-TOPAS parametric am- 
plifier (Light Conversion Ltd.). We use the idler from 
the TOPAS to get 800 /xJ pulses, continuously tunable 
between 1700 and 2200 nm. These pulses are focused by 
a 75 cm radius spherical silver mirror into a 2 mm contin- 
uous gas cell filled with argon (backing pressure around 
50 mbar). The pulse energy can be adjusted by rotating 
a 30 /im thick nitrocellulose pellicle, whose transmission 
varies from ~ 95% under normal incidence to ~ 50% un- 
der ^ 70° incident angle. The high harmonics are sent 
to an XUV flat field spectrometer consisting of a 250 /im 
slit, a 1200 mm~^ (variable groove spacing) grating (Hi- 
tachi), a set of dual microchannel plates associated to a 
phosphor screen (Hamamatsu), and a 12-bit cooled CCD 



The proper calibration of our XUV spectrometer is cru- 
cial to determine accurately the position of the Cooper 
minimum. The incidence angle of the grating is accu- 
rately set using a precision rotation mount. The dis- 
tances between the source and the grating and between 
the grating and MCP are measured. The central wave- 
length of the fundamental radiation is measured using a 
near infrared spectrometer. From all these parameters 
we can determine the theoretical position of the different 
harmonics on the MCP, assuming a certain order go for 
the first harmonic. We compare these theoretical posi- 
tions to the measured positions and determine the value 
of go which provides good agreement. Since we measure 
many harmonics, we can check the fine agreement of the 
calculated and measured positions over a broad spectral 
range. We can finely adjust the distance value between 
the grating and MCP, which is subject to the biggest un- 
certainty in our experiment, to achieve perfect agreement 
and extract the pixel- wavelength conversion function. As 
a check of the calibration, we measure the transmission 
of the aluminum foil in our experiment. We compare 
them to the theoretical transmission of a 200 nm thick 
filter with a 5 nm layer of alumine on each side [29| . The 
position of the abrupt cutoff in the filter transmission is 
in excellent agreement. 

A typical harmonic spectrum obtained at 1830 nm 
laser wavelength is shown in Fig. [5] The spectrum was 
normalized taking into account the pixel- wavelength con- 
version, the diffraction efficiency of the grating and the 
measured transmission of the aluminum filter. An ad- 
ditional advantage of using long wavelength appears in 
this spectrum: the fundamental photon energy is 0.68 
eV so that the spectrum is more dense than at 800 nm 
(to which corresponds a photon energy of 1.55 eV), mak- 
ing the determination of the position of the minimum 
more accurate. The dashed line in Fig. [5] is a gaussian 
smoothing of the spectrum, from which we extract the lo- 
cation of the Cooper minimum. We performed a series of 
identical measurements to evaluate the average position 
of the minimum and obtained E = 53.8 ± 0.7eV. 
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FIG. 2: High harmonic spectrum generated in argon by a 
1830 nm 8 x 10^^ W.cm~^ laser pulse. The spectrum is av- 
eraged over 25000 laser shots. The dashed line is a gaussian 
smoothing of the spectrum. The arrow indicates the position 
of the Cooper minimum E = 53.8 eV. 



B. Systematic study 

In order to check the robustness of the position of 
the Cooper minimum against experimental conditions, 
we have performed a systematic study of HHG in argon, 
varying both the laser field parameters, which control the 
single atom response, and the macroscopic parameters 
(gas pressure, beam focusing conditions), which influence 
the finally detected HHG yield. 



1. Laser field parameters 

An important difference between HHG and photoion- 
ization experiments is the presence of a strong laser field, 
which might modify the position of the Cooper minimum. 
We thus conducted a study as a function of the laser field 
parameters (intensity and wavelength). 

Figure El^a) shows a few harmonic spectra produced at 
different laser intensities (controlled by rotating the ni- 
trocellulose pellicle). As the intensity increases the cut- 
off is shifted from 60 eV to 90 eV. The Cooper mini- 
mum is clearly visible on all spectra. Figure shows 
the measured position of the minimum as a function of 
the laser intensity. We do not observe any systematic 
variation of the minimum location as the intensity is in- 
creased twofold even if the ponderomotive energy changes 
by more than 10 eV. The positions at lower intensities 
are slightly below average but the signal being smaller 
the error bars are more important. We also varied the 
central wavelength of the laser pulses delivered by the 
HE-TOPAS between 1800 and 1980 nm (Fig. ^c)). We 
do not observe any significant shift of the Cooper min- 
imum, which stays around 53.8 eV. These two observa- 
tions constitute strong indications of the lack of infiuence 
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FIG. 3: (a) High harmonic spectra generated in argon by a 
1830 nm laser pulse at intensities between 5 x 10^'^ W.cm"'^ 
(bottom) and 8 x 10^^ W.cm~^ (top). No aluminum filter was 
used for these spectra. The spectra are arbitrarily shifted ver- 
tically with respect to each other. Each spectrum is averaged 
over 25000 laser shots, (b) Position of the Cooper minimum in 
Argon as a function of the laser intensity at 1830 nm and (c) 
as a function of the laser wavelength. The horizontal dashed 
line marks the minimum location value E = 53.8 eV deter- 
mined previously. 



of the laser field on the recombination process. 



2. Phase matching 

While the Cooper minimum is a characteristic of high 
order harmonic emission from a single atom, our exper- 
iment measures the outcome of a macroscopic process. 
It is thus important to evaluate the possible influence of 
propagation effects in the measured spectra [30] . To that 
purpose, we varied several parameters which affect the 
phase matching conditions and the macroscopic buildup 
of the harmonic signal. First, we varied the backing pres- 
sure of the gas cell between 12 mbar and 160 mbar and 
did not observe any change in the position of the Cooper 
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FIG. 4: Position of the Cooper minimum in Argon at 1830 nm 
as a function of the backing pressure in mbar (a) , the laser 
beam aperture (b) and longitudinal position of the focus (c). 
The horizontal dashed line marks the minimum location value 
E — 53.8 eV determined previously. 



III. PHOTOIONIZATION AND RADIATIVE 
RECOMBINATION 

Our experimental study has shown that the Cooper 
minimum in HHG is located at 53.8 eV. This position 
is different from that measured in photoionization spec- 
tra (between 48 and 49 eV). Even though at first sight 
photoionization and recombination appear as strictly re- 
verse processes, which would lead to a simple conjugation 
relation between the associated transition dipole, the ex- 
perimental observations show a systematic shift. In this 
section, we perform a detailed theoretical analysis of the 
link between photoionization and radiative recombina- 
tion processes in order to explain this difference. 

Our calculations rely on the Single Active Electron 
(SAE) approximation [32l | where the description of elec- 
tron dynamics is restricted to that of a single Ar valence 
electron. The interaction of this electron with the nucleus 
and other electrons is reproduced, in the framework of a 
mean- field theory, by the model potential (33| : 



-Br 



+ {17- A) e 



-Cr 



(1) 



with A=5A, 5=1, C=3.682. This model potential ful- 
fills the correct asymptotic condition, i.e. V{r — )■ oo) — 
— 1/r, and the expected behavior at the origin, V{r — > 
0) = — 18/r. In [33], the parameters A, B and C were 
adjusted so as to provide the correct /p(Ar); we have 
further checked that V{r) yields accurate values for the 
energies of the Ar(?i^) excited states with n > A. 

The diagonalization of the hydrogen-like Hamiltonian 
H ~ + V in a, basis of even-tempered Slater- type 
orbitals (STO) yields bound and discretized continuum 
states in the spherical form ^„,;,m(r) — Rn^i{r)Yi^{V,) 
where YJ"* are spherical harmonics of given (I, m) symme- 
try and fir = {Or, (fir) is the angular f direction. The scat- 
tering continuum states ^I'k, normalized on the wavevec- 
tor scale k, are developed on the spherical state basis 



minimum (Fig. SKa)). Second, we changed the laser 
beam aperture between 11 and 17 mm. This resulted in a 
modification of the harmonic cutoff but no change in the 
position of the Cooper minimum (Fig. lUJb)). Finally, we 
varied the longitudinal position of the laser focus with 
respect to the center of the gas cell (Fig. IH^c)). This also 
led to a change of the cutoff position but the value of 
the minimum location remained the same. We note that 
we did not observe any signature of the presence of long 
trajectories in our experiment while scanning the focus 
position: the harmonics always appear as spatially col- 
limated and spectrally narrow, which indicates that the 
short trajectories would always be favored in our gener- 
ating conditions. 

These observations are different from what has recently 
been reported by Farrell et al. in a study of the Cooper 
minimum in Argon using 800 nm pulses [311] . 



oo I 
1=0 m=~l 

(-i)'"i^™(r!r)rr'"(f^k) 

(2) 

where f2k = {9ii_,ipk) corresponds to the angular k di- 
rection. Since the diagonalization of H in the STO un- 
derlying basis yields a coarse-grained representation of 
the continuum [s^ , [35j , we alternatively obtain the ra- 
dial parts Rk,i{r), and the phase shifts 6k,i (including 
both Coulonibic and short range components), through 
direct integration of the radial Schrodinger equation, us- 
ing the Numerov algorithm. It is important to note that 
eq.(l2]) explicitly includes the influence of the ionic core 
on the ejected/recombining electron, beyond the plane 
wave approximation which does not allow to faithfully 
describe neither photoionization nor HHG [1^ [s^ . 
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A. Photoionization cross section 



The photoionization rate induced by absorption of a 
radiation of amphtude Fq a nd p olarization axis n is given 
by the Fermi's golden rule |33 |: 



w;(k,n) = -fc|(*k| 2? 1*3,1,0) 



(3) 



where V = —FqU-t is the electric dipole moment op- 
erator. The photoionization differential cross section is 
given by: 

^!^^l![!fc^|(vl,k|„.r|vl,3^i,„)|2 (4) 

w being the frequency of the ionizing photon, related to 
the electronic momentum by the relation cu — Ip + k'^ /2. 

We aim to compare our calculations to photoioniza- 
tion measurements of Marr & West [2lj and Samson et 
al [22']. In these experiments, the measured quantity is 
the total photoionization cross section, which includes all 
possible relative orientation between the electronic mo- 
mentum and light polarization axis: 



a{k) = 



(5) 



Because of the selection rules for electric dipole tran- 
sitions, the calculation of (5'k| n • r |4'3,i.o) is reduced to 
the components of *k whose angular momenta are I — 
(p — > s transition) and I — 2 (p ^ d transition). The 
angular part of this calculation can be done analytically, 
and the total photoionization cross section can finally be 
expressed as: 



24%) 



(6) 



where c is the speed of light, and Ik,i are the radial inte- 
grals Ik.i = f Rk,i{r)R3,i{r)r^dr. These radial integrals 
are plotted in Fig. [5l //c,o is positive whatever k values, 
and slowly decreases as the electron energy increases. 
The p ^ d integral exhibits stronger modulations and 
a sign change for an electron kinetic energy of 30.6 eV 
(corresponding to a photon energy of 46.3 eV). This sign 
change is responsible for the minimum observed in the 
calculated total photoionization cross section, as seen in 
Fig. E 

Our calculation shows a clear minimum in the pho- 
toionization cross section for a photon energy of 50.1 eV, 
slightly above the measured value (between 48 an 49 eV) . 
The agreement with the experiment is fairly good, except 
for low photon energies (below 30 eV), as usual in single 
active electron calculations. In order to perform accu- 
rate calculations in this energy range, one has indeed to 
take into account multielectronic effects and transitions 
to excited states of the neutral 1371 . 
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FIG. 5: Radial integrals Ik,i associated with p ^ s (dashed 
line) and p ^ d (plain line) transitions. 
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FIG. 6: Calculation of the total photoionization cross section 
(plain line) compared to measurements from [SJ] (dotted line) 
and (dashed line). 



B. Recombination dipole in HHG 

1. Linear laser field 

An important specificity of high order harmonic gener- 
ation compared to photoionization is the selection of spe- 
cific values of k and n. First, tunnel ionization selects the 
quantization axis of the atomic orbital (z) parallel to the 
polarization axis of the IR generating field, as recently 
confirmed experimentally [38-40]. The IR field drives 
the trajectory of the recombining electron along the same 
axis. The electronic momentum k thus has to be taken 
parallel to the quantization axis of the orbital. Further- 
more, by symmetry consideration, the emitted XUV field 
has to be parallel to the common axis of the laser polar- 
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FIG. 7: Calculation of jdrecl'^ (plain line) compared to total 
photoionization cross section (dashed line). 

ization and quantization. The calculation is thus reduced 
to the case of k |j n || z. 

Within this frame, the recombination matrix clement 
is equal to: 

|drcc|' = |(*k|n.r|*3,i,o)|' 

= J^l^k^oe^'"-" - 2/fe,2e'''='1' (7) 

The evolution of |drccP as a function of the emitted 
photon energy is plotted in Fig. [T] The dipole moment 
shows a Cooper minimum which is more contrasted than 
that observed in the total photoionization cross section. 
This is due to the coherent nature of the differential cross 
section: the differential calculation leads to a coherent 
sum of the p — > s and p ^ d transitions, which depends 
on the values of the scattering phases Sk.o and Sk.2- In 
the case of angular integrated calculation, this phase in- 
fluence vanishes leading to an incoherent sum of the two 
contributions. 

The Cooper minimum in the recombination dipole mo- 
ment is located at 51.6 eV, while it is at 50.1 eV in the 
total photoionization cross section. The selection of the 
quantization axis and light polarization thus lead to a sig- 
nificant shift of the Cooper minimum. In the following, 
we confirm the importance of this selection by consider- 
ing the case of an elliptical laser field. 

2. Elliptical laser field 

The recoUision direction of the electron in HHG can 
be manipulated using an elliptically polarized laser field. 
In that case, the electron trajectory between tunnel ion- 
ization and recombination is two-dimensional, the recol- 
lision axis is not parallel to the quantization axis and the 
polarization direction of the harmonics is not that of the 
laser field ^J. 

We have performed classical calculations to study the 
influence of ellipticity on the Cooper minimum, in the 



framework of the so-called " Simpleman" model [l| that 
neglects the influence of the ionic potential on the elec- 
tron dynamics between tunnel ionization and recombi- 
nation. The electron trajectory is supposed to be ion- 
ized at time ti through tunneling along the instanta- 
neous field direction which sets the quantization axis; 
this accordingly yields the initial conditions r{ti) = 
and W||(ti) = 0. We further require the trajectory to be 
closed, i.e. r{tr) = where tr is the recombination time. 
This imposes a non-zero perpendicular velocity at time of 
ionization, v±{ti) ^ 0, consistently with the lateral con- 
finement of the electronic wavepacket during tunneling. 
The recoUision angle 9k, which is defined with respect 
to the quantization axis and determines the k direction 
at time of recombination '^\, is then collected for each 
electron trajectory. This finally allows us to compute the 
parallel d|| and orthogonal d±^ components of the recom- 
bination dipole. 

In Fig. [5Ja), we represent the value of |(ircc(fc, ^k)P = 
\d\\ + d±\'^ as a function of the generating field elliptic- 
ity e, for A = 1900 nm and / =lxlO"W.cm"2. As e 
increases the overall shape of the minimum is modified: 
it becomes broader on the high energy side, with an ap- 
proximate slope of 0.15 eV by percent of ellipticity. This 
behavior stems from increasing contributions of 0k 7^ 
recombination angles to the HHG signal as e increases, 
given that the minimum of |(ircc(fc, 0k)P shifts to higher 
energy as 0^ increases. 

We have performed measurements to confirm the 
broadening of the Cooper minimum. We have added 
a quarter waveplate in the laser beam. The ellipticity 
was controlled by rotating the axis of the waveplate with 
respect to the laser polarization. As e increases the har- 
monic signal falls rapidly, preventing measurements with 
ellipticities larger than 20%. Figure [5Ub) shows the nor- 
malized harmonic spectra as a function of the laser ellip- 
ticity, at an intensity of 1 x lO^'^W.cm"^ and a laser wave- 
length A=1900 nm. The white dotted lines are identical 
in the two panels, which enables a direct comparison. 
As predicted by our simulations, the Cooper minimum 
broadens on the high energy side as we increase the e, 
with an average slope of approximately 0.2 eV by per- 
cent of ellipticity. These results confirm the important 
role of the recoUision direction in the accurate determi- 
nation of the shape of the harmonic spectrum even in a 
simple atomic system like Ar. 

Even though the calculated dipole moments and mea- 
sured harmonic spectra show very similar shapes, there 
is still a systematic 2.2 eV shift between the positions 
of the measured and calculated Cooper minima. This 
is due to the fact that the harmonic spectrum S{io) is 
not only determined by the recombination dipole moment 
|c?roc(fc, ^k)P but also by the number of recombining elec- 
trons at each energy, i.e. by the shape of the recoUiding 
electron wavepacket. In the following, we use a CTMC 
approach to take this effect into account and compare it 
to the result of our experiment. 
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FIG. 8: (a) Calculation of ]dioc(fe, Sk)P as a function of the 
photon energy and the ellipticity of the driving field (A=1900 
nm, .cm ). (b) Measured harmonic spectra in 

the same conditions. The dotted white lines are the same in 
the two panels. 



IV. COMPLETE THEORETICAL 
DESCRIPTION: CTMC-QUEST 



target system effectively considered. Our semi-classical 
description aims at avoiding such a requirement of ap- 
proximated or additional computations; we therefore use 
CTMC simulations to build the returning wavepacket 
in terms of electron trajectories which fulfill preselected 
rescattering conditions. Once the wavepacket is defined, 
quantum rescattering theory is used to compute the re- 
combination probability. 

A. CTMC and initial phase-space distribution 

The CTMC method has originally been developed to 
describe non-adiabatic processes in atomic collisions [1^ . 
It has been successfully applied to shed light on sub- 
tle ionization mechanisms (45l . |46| that cannot be rep- 
resented unambiguously by purely quantum mechanical 
treatments. CTMC calculations have also been per- 
formed to describe laser-matter interactions see e.g. [43- 
[igj . As in atomic collisions, the classical assumption in- 
herent in the CTMC method is not prohibitive in the de- 
scription of field-induced processes; for instance, CTMC 
trajectories leading to HHG have been shown to match 
their quantum counterparts obtained in the framework of 
an hydrodynamical Bohmian description of the electron 
fiow [50]. 

The CTMC approach employs a A/'-point discrete rep- 
resentation of the phase-space distribution g{r, p, t) in 
terms of independent electron trajectories {rj{t),pj{t)} 

1 ^ 

i=i 

where p is the canonical momentum conjugate to r. The 
temporal evolution of the g distribution is governed by 
the Liouville equation, which is the classical analogue to 
the time-dependent Schrodinger equation 



In order to get a full description of the HHG process 
that includes the influence of the structure of the recol- 
liding electron wavepacket, we have developed a semi- 
classical theoretical simulation called CTMC-QUEST 
(Classical Trajectory Monte Carlo-QUantum Electron 
Scattering Theory). Our theory is based on a combi- 
nation of classical and quantum approaches in the mod- 
elling of the HHG process. It is largely inspired by the 
well known three-step model of laser-matter interactions 
[H , in which most of the processes subsequent to primary 
field-induced ionization are determined by the recoUision 
of the returning wavepacket on the ionic core. In this re- 
spect, it is noteworthy that a Quantitative Rescattering 
Theory (QRS) has recently been developed to success- 
fully describe, e.g., HHG and non-sequential double ion- 
ization processes based on the recollision picture [13, . 
In the QRS, the returning electron wavepacket is ob- 
tained by means of the strong field approximation (SEA) 
[il ] or time-dependent Schrodinger calculations of an eas- 
ily solvable model system with similar Ip than that of the 



dg{r,p,t) dg{r,p,t)dH dg{r,p,t) dH 



dt 



dr dp 



dp dr 



(9) 



where H = H + v ■ F(t) = + V{r) + r • F{t) is the 
total electronic Hamiltonian, including the laser-target 
interaction r • F{t) expressed in the length gauge. In 
our calculations, we use a simple sinusoidal laser field 
F{t) = i^osin(a;t) where t € [0,47r/a;]. This description 
reduces the interaction of a real shaped laser pulse with 
a macroscopic medium to a simplified scheme where the 
harmonic signal is generated by a single atom submitted 
to an effective laser field intensity. This simplification is 
justified by the insensitivity of the experimental results 
to macroscopic effects and to the laser intensity. 

Inserting eq.® into eq.® yields the Hamilton equa- 
tions which tailor the motion of the j*'^ classical trajec- 
tory among the set of J\f independent ones 



dt 
dPjjt) 
dt 



-Vr(F(r) 



m)Ut) 



(10) 
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The last equation emphasizes how the electron evolves 
under the combined action of the ionic and laser fields, 
beyond the widely used SFA that neglects V{r); this has 
important consequences on the description of radiation 
and electron emission processes close to the ionization 
threshold 

The most common way to define an initial distri- 
bution is to use a microcanonical distribution whose 
energy is —Ip. In Fig. |9l we compare the quan- 
tum radial probability density 47rr^i?| ^^(r) of the ini- 
tial 3p Ar state to its microcanonical counterpart 
PM{~Ip;r t = 0) = / / Qm{~Ip;r,p,t ^ 0)dpdnr = 
167T^r^Cy/2{~Ip + V{r)) where C is a normalization con- 
stant. Besides the fact that pM{r,t = 0) does not repro- 
duce the nodal structure of the quantum i?§ ^ ('') j it re- 
stricts the electron density to an inner region close to the 
nucleus (r < 2.5) where ionization processes are classi- 
cally less probable. Such distribution can thus not prop- 
erly describe the ionization. 

Better quantitative results can be obtained in 
the framework of the CTMC approach provided 
that one employs an improved initial distribution 
g(r,p,t = 0), beyond the usual microcanonical descrip- 
tion gm{—Ip',J^,P,t = 0) which affects —Ip to the energy 
of all Af trajectories [5J| . 

Improving the initial condition is based on the obser- 
vation that a trajectory with energy Ej{t — 0) — pj{t ~ 
0)+V{rj{t — 0)) ^ —Ip can participate to the description 
of the initial state since the negative part of the classical 
energy scale is not quantized. In other words, one can 
construct an improved initial phase-space distribution as 
a functional integral over microcanonical ones 



^?(r,p,t = 0)= / c{E)QM{E;r,p,t^O)dE (11) 

JE- 

where the bounds E^ and E^ are given by the parti- 
tion of the classical phase space into adjacent and non- 
overlapping energy bins [E~ ,E^[ for given I. The parti- 
tion is made following the method explained in [5J|; for 
Ar(3p) this leads to E' = -1.885 and E+ = -0.150 for 
a classical momentum |L| = |r x p| enclosed in [1,2[. In 
practice, the improved g{r, p,t — 0) consists of a discrete 
representation of ea.(|lip in terms of 10 microcanonical 
distributions; the (discretized) coefficients c^;^ are ob- 
tained by fitting the quantum radial probability density 
to p{r, t = 0) = Ylm=i CE^PM{Em;r, t = 0) with the ad- 
ditional condition that < E > p— —Ip. 

The final p{r, t — Q) \s displayed in Fig. |9l while the 
nodal structure of the quantum density still remains be- 
yond the scope of the improved classical description, as 
expected, the newly defined p{r,t = 0) nicely matches 
the outer region of i?| ^ (r) from which classical electrons 
preferentially escape subject to the laser field. Finally, 
our improved calculations employ A/" = 20 x 10^ trajec- 
tories in order to minimize statistical uncertainties. 



T p— -1 \ ^ I ^ \ ^ \ r 



0.6 — 




1 2 3 4 5 

r (a.u.) 

FIG. 9: Comparison between the initial density p(r, t = Q) as 
a linear combination of 10-microcanonical distribution (plain 
line), pM{—Ip',r,t = 0) the one-microcanonical distribution 
(dashed line) and the quantum radial electronic density of 
the fundamental Ar 3p state (dotted line). 

B. CTMC and returning electrons 

We now address the definition of the rescattering elec- 
tron wavepacket. We define a rescattering sphere, cen- 
tered on the target nucleus, of radius i?rcc of the order 
of the extension of the fundamental wavefunction. We 
consider that an electron is rescattering, and will gener- 
ate an harmonic photon through recombination, if after 
being ionized and leaving the sphere, it comes back into 
it. The Fig. [TU] shows a typical short electron trajectory 
which contributes to the harmonic signal. The recombi- 
nation time is taken as the time the electron enters back 
into the sphere. We record the energy Ej and the di- 
rection k of the wavevector at this instant. Within the 
CTMC statistics, we are thus finally able to define the 
density of the returning wavepacket at time t 

A/-.et(t) 

gretiEXt) = j^ SiE,{t)-E)SCk,it)-k) (12) 

i—l 

where Mret{t) is the number of electron trajectories 
that fulfill at time t the rescattering criterion mentioned 
above. 

Figure [TT] shows a typical result for / — 10^"* W.cm"^ 
and A — 1830 nm. By monitoring the energy of the 
returning electron with respect to the recombination time 
(Fig. [Ula)) we can identify the short (plain line) and 
long (dotted line) trajectories. This is an important asset 
of CTMC-QUEST: the separation of the different classes 
of trajectories (short, long, or even multiple return ones) 
is natural. 

In Fig. [HJb), we plot the time- and angle-integrated 
density of returning electrons for long and short trajec- 
tories 

Pret[E)= j Qret{EXt)dUt (13) 
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FIG. 10: Typical short electron trajectory that contributes 
to the harmonic signal. The electron after being ionized and 
leaving the sphere comes back into it. The inset shows the 
recombinaison angle 9k when the electron is recombining. 



The density of returning electrons pret{E) increases 
with energy for fong trajectories whereas it decreases 
for short ones. In other words, pret{E) decreases as the 
trajectory length increases. Two opposite effects are at 
play to determine the evolution of pret {E) with trajectory 
length. First, longer trajectories are born at earlier times 
when the laser field is more intense and thus benefit from 
a stronger ionization rate. However longer trajectories 
also show a stronger lateral spread during propagation 
in the continuum, which reduces the number of recom- 
bining electrons. Our study shows that the latter effect 
is dominant. It is consistent with the fact that we do not 
observed long trajectories in our experimental study. 

In Fig. [THwe present pret{E) for several values of the 
photorecombination sphere radius -Rrcc- This figure en- 
ables us to check the convergence of our calculation. In 
the plateau region, there is no significant difference be- 
tween the shapes of the wavepackets issued from the 3 
calculations, which means that the convergence is good. 
In the following, we will use a photorecombination sphere 
of -Rrcc = 5 a.u., which roughly corresponds to the exten- 
sion of the fundamental Ar wavefunction (see Fig. 

Our results can be compared to those obtained by 
Levesque et al. [11] who used experimental spectra and 
dipole moments calculated within the plane wave approx- 
imation (PWA) to extract the amplitude of the recollid- 
ing electron wavepacket. The lower wavelength of the 
driving field (800 nm) employed by Levesque et al. can- 
not explain the sharp disagreement between their and our 
results: while the shape of our wavepacket falls down by 
(only) a factor of ~ 5 in a 100 eV energy interval, several 
orders of magnitude appear between the low and high 
energy signals in [sl]. We checked that this discrepancy 
is due to the use of PWA by dividing our experimental 
HHG spectrum of Fig. [Hby the square of the dipole mo- 
ment computed in the PWA and recovering the drastic 
decrease of pret{E) as E increases. This confirms the 




FIG. 11: (a) : Total energy of the returning electron as func- 
tion of the photorecombination time. This representation per- 
mits a direct and intuitive separation of the different trajec- 
tories. In this figure, we only represent the long and short 
trajectories and omit the multiple returnings. (b) : Repre- 
sentation of the density of returning electrons for each type 
of trajectory. 



importance of taking into account the influence of the 
ionic potential in all the three steps of the laser-matter 
interaction. 

We can also compare our results to the Quantitative 
Rescattering Theory (QRS) calculations from Le et al. 
(see Fig. 2 of In the QRS framework, the ampli- 

tude of the electron wavepacket slightly depends on the 
method employed to compute the HHG spectrum. Le 
et al. display a flat-shaped wavepacket amplitude, as 
a function of E, when TDSE calculations are used to 
compute the spectrum before dividing it by the dipole 
moment. Significant deviations are obtained at low E 
when the SFA is employed. We have performed (but do 
not show for sake of conciseness) similar calculations to 
those of Le et al. who employed a 800 nm driving pulse. 
This has allowed us to ascertain that at 800 nm the shape 
of pret{E) is indeed flat from low to high E when inte- 
grating for all trajectories and multiple returns as it is 
the case in TDSE calculations. 
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FIG. 12: Comparison of the convergence for different radius 
size Tirec = 3 (plain line), 5 (dashed line) and 7 a.u. (dotted 
line) of the photorecombination sphere for the short trajecto- 
ries. 



C. QUEST 

The photorecombination probabihty PpB{E,'k,t) is 
quantum mechanically obtained by using the microre- 
versibility principle and the Fermi's golden rule : 

= ^^l(*^.,k|n-r|*3,i,o)P (14) 

where ^ is the stationary scattering wavefunction for 

a total energy E and a wave vector direction k. In the 
case of photoionization, the electron wavevector k' is re- 
lated to the kinetic energy at infinity through E = k'^/2. 
Recombination occurs at short range where the total elec- 
tron energy reads E = + V{Rycc)- The analogy 
between photoionization and photorecombination is thus 
only respected if one sets k' = \plE. 

In the first part of the manuscript, we calculated pho- 
torecombination probabilities assuming that the electron 
trajectory was linear and parallel to the quantization 
axis. In the present statistical 3D modelling, this is no 
more the case: the electron escapes from r 7^ with 
a non-zero transverse velocity; its interaction with the 
ionic potential can further deviate its trajectory from a 
straight line. Therefore, it is necessary to take the re- 
combination angle into account, defined with respect 
to the z-axis, in the computation of the photorecombina- 
tion probability 

\FXE 

Vpn{EM = -^\lE,i..^^''^'^'' 

(3cos2 0fe-l) e'*-^M|2 
(15) 



FIG. 13: Comparison between the statistical distribution of 
recombination probabilities (scattered points which look like 
a plain line, see inset for details) and the harmonic spec- 
trum (dashed line) for short trajectories obtained by CTMC- 
QUEST. The position of the Cooper minimum is shifted to- 
wards high energy for the harmonic spectrum. This is a con- 
sequence of the shape of the density of returning electrons. 



D. CTMC-QUEST for HHG 

Once the density of returning electrons and photore- 
combination probabilities are determined, we can calcu- 
late the harmonic spectrum 5(0;, n,t) as: 

5(w,n,t)^ y" dEerec(EM,t)VpR(EM)8{E^l^-u:) 
= 77 E VpRmt)Mt))f{Em,uj) (16) 

i—l 

where f{Ei{t), w) = 1 if Ei{t) = uj — Ip and otherwise. 
In case of a laser field linearly polarized along the z-axis 
n — z . The time-integrated HHG spectrum is simply 
obtained through S{uj,h,T) — S{uj,h,t) where r is 
the pulse duration. 

1. Role of the returning electron wavepacket 

In Fig. [131 we display the statistical distribution of the 
recombination probabilities VpR{Ei{t),'ki{t)) associated 
to all the short trajectories that lead to HHG during 
the whole interaction. This distribution consists of a set 
of scattered points in the {E,'PpR{E,'k)) plane (see the 
inset in Fig. I13p . The distribution is very narrow since 
dki ~ for most of the trajectories so that its shape is 
similar to that of the dipole of Fig. [T] 

The time-integrated HHG spectrum S{oj, n, r) built ac- 
cording to eq. (1161) with only short trajectories is included 
in Fig. 1131 The position of the Cooper minimum in 
5(a;, n, r), located at £^ = 53.5 ± 0.5 eV, is different from 
the one obtained in the probability distribution. The 
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FIG. 14: Comparison between the theoretical harmonic spec- 
trum obtained either for long (dashed line) and short (plain 
line) trajectories. As expected looking to the shape of the 
returning wavepacket, the position of the Cooper minimum is 
different for these two types of trajectories. 
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FIG. 15: Comparison between an experimental (dashed line) 
and CTMC-QUEST (full line) spectra obtained for with a 
1830 nm laser field at 1 x lO^^'W.cm-^ 



latter one, found at E ^ 51.7 ± 0.5 eV, is close to the 
value obtained assuming 0^. — for all the trajectories. 
The difference of minimum locations is mainly due to 
the shape of the recoUiding electron wavepacket for short 
trajectories (Fig. [TIT b)): the decrease of pret{E) with in- 
creasing E shifts the minimum observed in the photore- 
combination probability distribution to higher energies. 
This result, obtained on a simple atomic system, demon- 
strates that the accurate study of orbital structure via 
HHG cannot be restricted to the study of the dipole mo- 
ment. The shape of the returning wavepacket has also to 
be accurately known. 

The importance of the shape of the recoUiding 
wavepacket is confirmed in Fig. [TJ] where we compare the 
harmonic spectra generated by the short and the long tra- 
jectories. Interestingly, the position of the Cooper mini- 
mum is sensitive to the considered set of trajectories: it 
stands at 51.7 eV for long trajectories and at 53.7 eV for 
short ones. This behavior is due to the opposite slopes 
of the electron wavepacket profiles (see Fig. [TlTb)). 

2. Comparison with experiments 

Figure [15] shows a comparison between the experimen- 
tal and theoretical HHG spectra for A = 1830 nm and 
/ = 10^"* W.cm^-^. The agreement is very satisfactory 
since the position of the theoretical Cooper minimum is 
53.7 eV which nicely matches the experimental one, 53.8 
±0.7 eV. The overall shape of the experimental spectrum 
is also very well described. 

In Fig. [ini we compare the results of our CTMC- 
QUEST calculations for three different laser intensities 
1 = 5, 7.5 and 10 x 10^^ W.crn"^. There is no clear vari- 
ation of the position of the Cooper minimum between 
the two highest intensities. However, at the lowest in- 
tensity the minimum is shifted to lower energies (52.7 
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FIG. 16: CTMC-QUEST spectra obtained with a 1830 nm 
laser field at 5 x 10^^ (plain fine), 7.5 x 10^^ (dashed line) and 
1 X 10^" W.cm"^ (dotted line). 

eV). This reflects the intensity dependence of the shape 
of the rescattering wavepacket. Our experiments are con- 
sistent with the lack of variation of the minimum between 
7.5 X 10^^ and 1.0 x 10^"* W.cm"-^. At lower intensities we 
observed a slight downshift of the minimum which could 
be the signature of the modification of the wavepacket 
described by the theoretical results. 



V. CONCLUSION 

In this paper we have studied the link between pho- 
toionization and high harmonic generation by focusing on 
the Cooper minimum in Argon. This structural feature 
is common to the two processes but while the position of 
the minimum is observed between 48 and 49 eV in total 
photoionization cross sections, we measured it at 53.8 eV 
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in HHG. By performing a systematic experimental study 
using 1800 nm laser pulse, we have concluded that the 
position of the Cooper minimum was independent on the 
gas pressure and focusing conditions and to a large extent 
to the laser intensity. 

The observed shift of the Cooper minimum is partly 
due to the fact that the recombination process involved 
in HHG is the inverse process of a particular case of pho- 
toionization, in which the quantization axis of the atomic 
orbital, electron ejection direction and XUV photon po- 
larization are parallel. We have checked the influence of 
these parameters by manipulating the electron trajectory 
using ellipticity, and have observed similar effects in re- 
combination dipole moments and experimental harmonic 
spectra. 

The dipole moment is not sufficient to determine accu- 
rately the harmonic spectrum: an additional shift of the 
Cooper minimum is due to the structure of the recom- 
bining electron wavepacket. 

In order to reproduce the experimental observations, 
we have developed a model based on the combination 
of CTMC and Quantum Electron Scattering techniques. 
CTMC-QUEST takes into account the role of the ionic 
potential on the 3D trajectories of the electrons in the 
continuum. The recombination dipole moment corre- 
sponding to each individual electron trajectory is quan- 



tum mechanically computed taking into account the rec- 
ollision angle with respect to the quantization axis. This 
procedure enables us to obtain a very satisfactory agree- 
ment between theoretical and experimental spectra. It 
further allows to unambiguously discriminate between 
short, long and multiple return trajectories contributing 
to the recollision wavepacket. 

Our results show that accurate high harmonic spec- 
troscopy can be performed using long wavelength lasers. 
The degree of precision reached by the experiment has 
allowed us to refine our theoretical modelling. We are 
now working on the extension of CTMC-QUEST to poly- 
atomic molecules in the perspective of achieving complete 
3D description of the generation process [l4[ . 
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